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ABSTRACT 

We revisit the so called "blazar sequence", which connects the observed bolometric lu- 
minosity to the shape of the spectral energy distribution (SED) of blazars. We propose that 
the power of the jet and the SED of its emission are linked to the two main parameters of the 
accretion process, namely the mass of the black hole and the accretion rate. We assume: i) that 
the jet kinetic power is proportional to the mass accretion rate; ii) that most of the jet dissipa- 
tion takes place at a distance proportional to the black hole mass; iii) that the broad line region 
exists only above a critical value of the disk luminosity, in Eddington units, and iv) that the 
radius of the broad line region scales as the square root of the ionising disk luminosity. These 
assumptions, motivated by existing observations or by reasonable theoretical considerations, 
are sufficient to uniquely determine the SED of all blazars. This framework accounts for the 
existence of "blue quasars", i.e. objects with broad emission lines but with SEDs resembling 
those of low luminosity high energy peaked BL Lac objects, as well as the existence of rel- 
atively low luminosity "red" quasars. Implications on the possible evolution of blazars are 
briefly discussed. This scenario can be tested quite easily once the AGILE and especially the 
GLAST satellite observations, coupled with information in the optical/X-ray band from Swift, 
will allow the knowledge of the entire SED of hundreds (and possibly thousands) blazars. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

Fossati et al. (1998) studied three complete sample of blazars: the 
Einstein Slew survey (Elvis et al. 1992), the 1-Jy samples of BL 
Lacs (Kiihr et al. 1981), and the flat-spectrum radio-loud quasars 
(FSRQs) extracted by Padovani & Urry (1992) from the 2-Jy sam- 
ple of Wall & Peacock (1985). The total number of studied blazars 
was 126, and 33 of these were detected in 7-rays by the EGRET 
instrument onboard the Compton Gamma Ray Obser-vatory. These 
blazars were divided into radio luminosity bins, and the luminos- 
ity in selected bands was averaged to form the SED representative 
of the blazars in each bin. It turned out that the division into ra- 
dio luminosity bins well matched the division into bins of bolo- 
metric luminosity, and that all spectra could be described by two 
broad peaks, in a vF v representation, the first at mm/soft X-rays 
frequencies, the second in the MeV-GeV band. More importantly, 
a sequence appeared: blazars with greater bolometric luminosity 
had "redder" SEDs (i.e. smaller peak frequencies: LBL in the ter- 
minology of Padovani & Giommi 1995), and the high energy peak 
was more prominent. Blazars of lower bolometric luminosity had 
instead a "blue" SED (HBL in the terminology of Padovani & 



Giommi 1995) with the two peaks having approximately the same 
luminosity. 

This spectral sequence was interpreted by Ghisellini et al. 
(1998) as due to the larger radiative cooling suffered by the emit- 
ting electrons of blazars of larger power, where the radiation energy 
density seen in their comoving frames received a large contribution 
by photons produced outside the jet. As a model, they adopted a 
simple leptonic one-zone synchrotron inverse Compton model. In 
this scenario a stronger cooling resulted in a particle energy distri- 
bution with a break at lower energies (producing smaller peak fre- 
quencies) with more power emitted through the inverse Compton 
process (hence the dominance of the high energy peak). This pic- 
ture was later confirmed (Ghisellini, Celotti & Costamante, 2002; 
Celotti & Ghisellini 2008) when the new generations of Cherenkov 
telescopes allowed the detection of an increasing number of low 
power BL LacsQ. Through the modelling, the above studies found 
that there is a good correlation between the energy 7 pC akTieC 2 
of the electrons emitting at the peaks of the SED and the total 
(i.e. magnetic plus radiative) energy density as seen in the comov- 
ing frame of the blazar. Therefore the so called "blazar sequence" 
comes in two kinds: i) a purely phenomenological sequence, relat- 
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ing the SED shape with the bolometric observed luminosity, and ii) 
a more "theoretical" one, relating 7 pca k to the amount of radiative 
cooling. It is clear that the first kind is associated with observed 
properties: since it considers the brightest blazars, it is likely that it 
corresponds to the most aligned sources: for blazars that are (even 
slightly) misaligned, the observed luminosity of an intrinsically 
powerful source becomes smaller and the observed SED becomes 
"redder", contrary to the general sequence. Furthermore, the phe- 
nomenological blazar sequence had a rather incomplete SED cov- 
erage (only 33 out of 126 sources were observed in the 7-ray band, 
at that time). Thus, as it stands, the phenomenological blazar se- 
quence probably describes active states of the sources, not neces- 
sarily their averaged status. 

The second kind, instead, treats comoving quantities and is 
less permeable to these effects. Since any reasonable blazar lumi- 
nosity function predicts that the number of blazars increases for 
decreasing luminosities, one simple prediction of the "theoretical" 
scheme is that there should be more "blue" blazars than "red" ones. 

A critical review of the blazar sequence is presented by 
Padovani (2007), who pointed out three main tests that the blazar 
sequence should pass: 

(i) existence of an (anti-)correlation between the synchrotron 
peak frequency ^ s , pca k and the bolometric observed luminosity; 

(ii) non-existence of "blue" powerful objects; 

(iii) "blue" sources should be more numerous than "red" ones. 

Point (i) concerns observed properties: since red low luminos- 
ity blazars could be slightly misaligned sources, the existence of 
red blazar with small observed luminosities is not invalidating the 
blazar sequence, that even predicts them. The other two points are 
more important. In general, blazars with emission lines (FSRQs) 
have larger bolometric luminosities and "red" spectra. However, 
there are exceptions to this general rule. One example is RGB 
J1629+4008 (at a redshift 2 = 0.272), discussed in Padovani et al. 
(2002). This blazar has broad emission lines with equivalent width 
typical of FSRQs (~80 A), and its synchrotron spectrum peaks at 
~10 1B Hz, that corresponds to a "blue" SED. The luminosity of its 
broad emission lines is relatively modest (~7xl0 43 erg s _1 ); the 
modelling yields parameters consistent with the "theoretical blazar 
scheme" (i.e. its 7 pea k belongs to the correlation found in Ghis- 
ellini et al. 1998; see Fig. 9 of in Padovani et al. 2002). Another 
object showing (weak) broad emission lines and a blue SED is RX 
J1456.0+5048 (z = 0.479f]. As discussed in Maraschi et al. (2008) 
this source is very similar to RGB J1629+4008, and belongs as well 
to the 7 P eak(C^') correlation. 

Padovani et al. (2003) searched for other blue quasars using 
two large sample of blazars, the Deep X-ray Radio Blazar Sur- 
vey (DXRBs), and the ROSAT All-Sky Survey-Green Bank Sur- 
vey (RGB), for a total of about 500 sources. They found a rel- 
atively large number of possible candidates, but this finding was 
based on a rather poor characterisation of the SED, parametrised 
through the radio, optical and X-ray fluxes and the broad band 
spectral indices connecting the radio to optical (aRo), the opti- 
cal to X-ray (aox) and the radio to X-ray (orx) fluxes. Note 
that in the absence of multi-band optical data it is rather difficult 
to disentangle the beamed non-thermal to the accretion disk blue 
bump continua. For this reason these results were interesting, but 
uncertain and therefore not conclusive. A small fraction of these 
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blue quasar candidates, observed at the VLA (Landt, Perlman & 
Padovani 2006), showed a rather modest core radio luminosity, at 
the boundary of the FR I and the FR II radio-galaxy division. Of 
these, 10 sources were observed by Chandra (Landt et al. 2008) 
and showed flat X-ray spectra (i.e. energy spectral index ax < 1), 
demonstrating that these blue quasars candidates have instead a red 
SED. Also the Sedentary survey (e.g., Giommi et al. 2005), tuned to 
find high synchrotron peak sources, has not detected flat spectrum 
radio quasars. 

Recently, two other high redshift and high power blazars 
were discovered, claimed to be blue quasars candidates: SDSS 
J08 1009.94+384757.0 (z = 3.95; Giommi et al. 2007) and IGR 
J22517+2218 (2 = 3.668; Bassani et al. 2007). Both are charac- 
terised by a large X-ray to optical ratio, that could be interpreted as 
due to a single synchrotron component peaking at X-ray frequen- 
cies. This possibility, mentioned in both papers, is disfavoured by 
the fact that the X-ray spectral slope, in IGR J22517+2218, is flat- 
ter than aox, suggesting that the X-ray flux belongs to the high 
energy peak. Instead, in SDSS J081009.94+384757.0, the optical 
flux can correspond to the thermal emission from the accretion disk, 
and if so, also this source fits very well in the blazar sequence (see 
the SED and the corresponding model in Fig. 3 of Maraschi et al. 
2008). 

The two papers mentioned above proposed also that if these 
sources are indeed red, then they would be the prototype of a new 
class of blazars with extreme properties, namely a very dominat- 
ing high energy peak. However, these two sources have a SED 
that is not unprecedented: the high redshift (z > 4) blazars (GB 
1428+4217, 2 = 4.72: Fabian et al. 1998; PMN J0525-3343, 
2 = 4.4: Fabian et al. 2001; RX J1028.6-0844, 2 = 4.276: Yuan 
et al. 2000: Q0906+6930, 2 = 5.47: Romani et al. 2004) have a 
very similar SED, with a similar X-ray to optical flux ratios (see 
also RBS 315 at z = 2.69; Tavecchio et al. 2007). Therefore all 
these sources are not powerful blue quasars, but the extreme man- 
ifestation of the blazar sequence (high energy peak increasingly 
dominating increasing the observed luminosity; see Maraschi et al. 
2008). Far from disproving the blazar sequence, they fully confirm 
it. 

Concerning the third test that the blazar sequence should 
pass (blue BL Lacs should be more numerous than red blazars), 
Padovani et al. (2007, see also the review of Padovani 2007) pointed 
out that the ratio of blue/red counts of BL Lacs in deeper (in flux 
limit) radio and X-ray samples of blazars disagrees with the pre- 
dictions made by the Fossati et al. (1997) assuming the blazar se- 
quence. The disagreement is moderate for the X-ray surveys, and 
is more severe for radio surveys down to the 50 mJy flux level (in 
the sense that Padovani et al. 2007 find a ratio blue/red BL Lacs a 
factor ~3 smaller than predicted by Fossati et al. 1997). 

The possible solution to this problem offered by the blazar se- 
quence is the following: if red blazars are intrinsically more pow- 
erful than blue BL Lacs, they can enter a flux limited sample even 
if the jet is (slightly) misaligned, while blue BL Lacs cannot. This 
selection effect could make red blazars to be over-represented in 
the sample. This is admittedly no more than an educated guess, and 
it should be supported by detailed simulations: Fossati et al. (1997) 
did not consider slightly misaligned blazars. As discussed below, 
the blazar sequence we are proposing in this paper can offer an al- 
ternative (or additional) explanation to the above problem, in terms 
of BL Lacs whose black hole has a relatively small mass. 

In this paper we try to improve and to extend the blazar se- 
quence, taking into account that 

i) the phenomenological sequence was based on bright sam- 
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pies, and on a few high energy detections, not on the "average" 
states of typical sources; ii) we now know many more TeV BL Lac 
objects; iii) we made progress in calculating the power of the jet 
of blazars; iv) we also made progress in estimating the mass of the 
central black hole in some blazars. 

The old phenomenological sequence was based on only one 
parameter: the bolometric apparent luminosity. We here explore the 
possibility to associate the SED of all blazars to the two fundamen- 
tal parameters of all AGNs: the mass of the black hole M and the 
accretion rate Mi B . To this aim we build upon the old "theoreti- 
cal" blazar sequence, driven by a few simple ideas and pieces of 
evidence: 

• There is a preferred region of the jet where most of the ob- 
served radiation is emitted (see e.e. Ghisellini & Madau 1996). The 
location of this region, in the jet, should scale as the mass of the 
black hole. 

• Broad emission lines come from a region whose distance from 
the black hole scales (approximately, and with some scatter) with 
the square root of the disk luminosity. 

• The broad emission line region exists only if the disk luminos- 
ity is above a critical value (in Eddington units). 

• The power of the jet (Poynting flux plus kinetic) scales as the 
mass accretion rate. 

We will examine and justify these points in the next section, 
here we stress only that these main hypotheses will suffice to com- 
pletely describe the SED produced by the jets of blazars. Therefore 
we will construct a two-parameter sequence, for which the SED 
is dependent upon M- m and M. Since the disk luminosity Ldi s k is 
univocally determined given Afj n and M, the two parameters can 
equivalently be Ldisk and M. 



2 ASSUMPTIONS 

In this section we discuss the main assumptions of our model. 



2.1 Black hole masses 

It is generally believed that the black hole mass of a radio-loud 
AGN is on average larger than the typical mass of a radio-quiet 
object. However, the exact range of black hole masses of radio- 
loud AGNs is still a matter of debate. 

Lacy et al. (2001) proposed that the radio-loudness is a func- 
tion of mass, with most radio-loud sources having M >a few 
xlO 8 M e , while D'Elia et al (2003), with a sample of strong 
lined blazars, found a range of masses extending significantly to 
the lower end, down to a few x 10 7 Mq (see Metcalf & Maglioc- 
chetti 2006 for a recent discussion). In Fig. Q] we show the distri- 
butions of masses of FSRQs in the D'Elia et al. (2003) sample and 
of BL Lac objects (Woo et al. 2005; Wagner 2008), and compare 
them with the range of black hole masses estimated by Ghisellini & 
Celotti (2001) for FR I and FR II radio-galaxies. A caveat is in or- 
der: the correlations used to calculate the mass (i.e. the mass-bulge 
optical luminosity, or the size of the broad line region and the line 
width, or the a — M relation) have a large scatter, and therefore the 
shown distribution are only approximate. We indicate in the figure 
the corresponding averages of log M and their la dispersion. We 
conclude that most black hole masses are in the range 3 x 10 7 - 
3 x 10 9 Mq, and use this range in the following. 




7 8 9 10 

Log M/M 



Figure 1. Distribution of black hole masses calculated in Ghisellini & 
Celotti (2001) for FR I and FR II radio-galaxies (top two panels) by D'Elia 
et al. (2003) for FSRQs and for BL Lacs by Woo et al. (2005) and Wagner 
(2008). In the bottom panel the shaded histogram is for TeV BL Lac only. 
The arrows indicate the average of the logarithm of the masses, as labelled. 



2.2 Accretion rates 

A natural upper limit is obviously the Eddington one. Another criti- 
cal accretion rate, probably connected with the change of accretion 
mode, is for accretion luminosities L c around a few x 10 _3 I/Bdd. 
Below this critical values the accretion flow can be advection dom- 
inated (Narayan & Yi 1994), or ADIOS (Blandford & Begelman, 
1999, 2004), becoming radiatively inefficient, hotter and geometri- 
cally thicker. In Ghisellini & Celotti (2001) we have proposed that 
the dividing line between FR I and FR II (in a radio luminosity 
vs host galaxy optical luminosity plot) corresponds to this critical 
value. There also are several evidences that BL Lac objects indeed 
have radiatively inefficient disks, whose thermal emission is never 
seen. The absence of a strong ionising luminosity produced by the 
disk explains in a natural way the absence (or the extreme weak- 
ness) of the emission lines in BL Lacs. Also, if this is the reason 
of the FR I vs FR II and BL Lac vs FSRQ behaviour, there are in- 
teresting consequences on the redshift evolution of these classes of 
sources, since it is conceivable that there is some evolution in the 
accretion rate (as in radio-quiet quasars), and therefore BL Lacs 
(low accretion rates) might be more local (no or negative evolu- 
tion), while FSRQs should be more numerous (and powerful) in the 
past (positive evolution), as proposed and discussed by Cavaliere 
& D'Elia (2002) and Bottcher & Dermer (2002), and further dis- 
cussed below. As a lower limit to the accretion rate we are guided, 
for instance, by M87, whose disk should emit at L ~ 10 _6 LEdd- 
Also with these very small accretion rates, radio-loud systems can 
produce powerful jets. 

To summarise, we will assume that for L c /L^dd < 
L/L^dd < 1 we have "standard" accretion disks which orig- 
inates the jets in FSRQs (and FR II radio-galaxies), while, for 
10 -6 < L/Lvdd < L c /L-Edd we have BL Lac objects (and FR I 
radio-galaxies). The exact value of this subdivision (i.e. Lc/^Edd) 
is not very important, as long as the subdivision exists. 
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2.3 Seed photons from the Broad Line Region 

For our purposes, one important parameter is the radiation energy 
density of the broad emission line photons. We therefore need to 
estimate the radius of the broad line region and its total luminosity. 
We assume that the latter is 10% of the accretion disk luminosity (if 
this is larger than L c , see above). For the radius of the BLR there 
are, in the literature, several proposals: 

• according to Kaspi et al. (2005), the relationship between the 
radius of the BLR -Rblr. and the ionising disk luminosity is: 

0.69±0.05 



-Rblr 
10 It days 



= (2.23 ±0.21) 



AL A (5100A) 



10 44 erg s- 



(1) 



where La(5100A) is the monochromatic luminosity calculated at 
5100 A. 

• Bentz et al. (2006) pointed out a source of uncertainty in this 
relation associated with the flux of the host galaxy, contributing 
more at low AGN luminosity. Considering then a sample of AGN 
observed by HST, they derived: 

log -Rblr = (-21.2±1.7)+(0.518±0.04) log [AL A (5100A)] (2) 

with Rblr in light days and AL A (5100A) in erg s _1 . 

• More recently, Kaspi et al. (2007) considered the CIV line and 
the continuum at 1350 A, deriving: 

. l 0.52±0.04 



Rblr 
10 It days 



(0.17 ±0.04) 



AL A (1350A) 



10 43 erg s- 



(3) 



Note that the last two relations have consistent slopes, but inconsis- 
tent normalisations if \L\ at 1350 and 5100 A is similar. 

Given the above uncertainties, we chose to assume the sim- 
plest hypothesis, which is a BLR radius scaling with the square root 
of the disk luminosity. Also, for simplicity, we use the bolometric 
disk luminosity, by assuming: 

-Rblr. = 10 17 L d f 2 k45 cm (4) 

This implies that the energy density of the line photons (for an ob- 
server at rest with the black hole) is constant: 

idisk 



E/blr = 0.1 



4 *"#blr c 



2.65 x 10 2 erg cm 3 



(5) 



where we have also assumed a covering factor equal to 10%. Fol- 
lowing Ghisellini & Madau (1996), we then assume that the ob- 
server comoving with the jet emission region measures (7 BLR given 
by 

E/blr = ^ T 2 IW = 3.76 x 10" 2 r 2 ergem" 3 (6) 

The spectrum of this component is the sum of Doppler broadened 
lines and a continuum. The most prominent contribution comes 
from the Lyo? line, and this spectrum can be well approximated by 
a blackbody, with a peak (in vF v , in the comoving frame) around 
2 x 10 15 T Hz (see Tavecchio & Ghisellini 2008 for a detailed dis- 
cussion). 

2.4 Seed photons from the pc-scale dusty torus 

In the unification scenarios for Seyfert 1 and Seyfert 2 AGNs, one 
assume the presence of a torus, at pc-scale distances, that blocks 
the broad line photons to observers at large viewing angles, and ab- 
sorbs low energy X-ray photons, making the received X-ray spec- 
trum very hard. Crucial to explain the X-ray background (Setti & 



Woltjer 1989; Madau, Ghisellini & Fabian 1994), the dusty torus 
could be present also in jetted sources, but probably only in FR II 
radio-galaxies (and FSRQ), since, in FR I radio-galaxies, the very 
nucleus is not hidden in the optical (Chiaberge, Capetti & Celotti 
1999). The possible importance of IR radiation from the torus for 
the inverse Compton emission of blazars has been pointed out by 
Blazejowski et al. (2000) (see also Sikora et al. 2002, 2008). 

The torus reprocesses the absorbed disk radiation into the IR 
band. The typical temperature is around 150-200 K (Cleary et al. 
2007), as indicated by recent Spitzer observations. We approximate 
the result of Cleary et al. (2007) by assuming that the torus repro- 
cesses half of the disk radiation (corresponding to a opening angle 
of 60°). The typical distances of the torus, Rm, is assumed to scale 
as £ d f 2 k , yielding a constant temperature (oc L/R m ). From the re- 
sult of Cleary et al. (2007) we then set: 



Rib, = 2.5 x 10 18 L H { 2 



(7) 



The corresponding radiation energy density, as measured in the co- 
moving frame, is 

f/i'a = 3 x 10~ 4 r 2 ergcirT 3 (8) 

The spectrum of this component is assumed to be a blackbody, with 
a peak in the comoving frame (in a vF v plot) at 3 x 10 13 T Hz. 

2.5 Location of the jet dissipation region 

We assume that the dissipation radius Rdi ss scales as the 
Schwarzschild radius R s : -Rdi ss = aR s , because this scaling is 
appropriate both in the case of an accelerating jet and in the sce- 
nario of internal shocks (Sikora, Begelman & Rees 1994; Ghisellini 
1999; Spada et al. 2001) where we have 



Rdis 



r 2 R 



a(r)R s 



(9) 



where Ro is the initial separation of two consecutive shells, which 
can be approximated as a multiple of the Schwarzschild radius. We 
allow ourselves to have a different factor a for BL Lacs and for 
FSRQs, to mimic the possibility to have two different origin for 
the main dissipation in these two classes of objects: namely, in BL 
Lacs, we could have a standing shock at some radial distance from 
the black hole (along the lines of Sokolov, Marscher & McHardy 
2004). Even in this case we assume that this distance is proportional 
to the black hole mass, but not necessarily with the same constant 
as in FSRQs. 

The size of the emitting region is assumed to be a cylinder 
of cross sectional radius r = i/> Rdiss and height AR' = r, as 
measured in the comoving frame. Here ip is the opening angle of 
the jet. 

2.6 The jet power 

We will base our considerations upon the results presented in 
Celotti & Ghisellini (2008). In that work we found that the jet 
power Lj in blazars is dominated by the protons associated to the 
emitting electrons. We showed that the jet Poynting flux is smaller, 
and also that the component of the jet power possibly transported 
by electron positron pairs cannot be dynamically important (see 
also Sikora & Madejski 2000, Maraschi & Tavecchio 2003). In 
blazars the thermal component due to the accretion disk luminos- 
ity is often hidden by the Doppler boosted non-thermal continuum 
of the jet, but in FSRQs the presence of the broad emission lines 
allows to estimate the disk luminosity. We found that the power 
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of jets is comparable, and often larger than the luminosity emitted 
by the accretion disk in FSRQs. This is even more true in BL Lac 
objects, where the disk radiation is almost always invisible. 

Based on these considerations, we propose the following 
ansatz: the jet power is always linked with the accretion rate, 
namely we can write, both for BL Lacs and FSRQs: 



(10) 



Since Ldisk = ??Mi n c 2 , the jet power can always be written as 

Lj — — Ldisk (11) 
V 

If we set, in FSRQs, the jet efficiency factor equal to the disk 
accretion efficiency r\ (in the standard accretion mode, i.e. ?)j = 
■q ~ 0.1), we have Lj = L disk . If, as it seems the case (Celotti & 
Ghisellini 2008; Maraschi & Tavecchio 2003) the jet power exceeds 
the disk luminosity, then > r\. 

We assume that the accretion rate in BL Lac objects is below a 
critical value, M c , at which the accretion flow changes regime from 
"standard" disk accretion to ADAF-like or ADIOS-like regimes, 
for which L d i s k oc M m (e.g., Narayan, Garcia & McClintock 
1997). 

We must then relate M n with the disk luminosity and then 
relate the disk luminosity with the jet kinetic power. Since at M c 
the two accretion regimes corresponds to the same Ldisk — L c , we 
have 



V = 0.1; 



M in > M c or Ldisk ^ L c 



V = 0.1 -^; M in < M c or L dis k «S L c 
M c 



(12) 



where we set L c — r\M c c . 

The above assumptions allows to consider Ldisk (or Lj) as the 
other key parameter, beside the black hole mass, instead of Mm. 
The obvious advantage is that both L d i s k and Lj can be derived (or 
even directly observed, in the case of Ldisk of FSRQs) in a much 
easier way than Mi n . 

The above ansatz implies that the outflowing mass rate can be 
found through 



Lj = rM out c 2 -> M out = %Mn 



(13) 



2.7 Jet Poynting flux and kinetic power 

We assume that a fraction es of Lj is in the form of Poynting flux. 
This allows us to estimate the magnetic energy density as: 

e B Lj 



Lb = £BLj = 7n/> 2 i?dissLsr 2 c 



Tn/^Ls^ 



(14) 



We assume that a fraction e c is converted, at i?diss, in relativistic 
electrons. Their energy density, calculated in the comoving frame, 
is U e = m c c 2 J N(y)dy and is given by 



e c Lj 



W> 2 i?iis S r 2 C 



(15) 



In fast cooling (i.e. when all the injected electrons can cool in one 
light crossing time) all the power in electrons is converted into ra- 
diation. In slow cooling, instead, only a fraction e\ can be emitted. 
Therefore we have that the energy density of the radiation produced 
by the jet is 

e e £ALj 



U T = CaU c = 



(16) 



In our scenario e c and es are free parameters, while ea can be de- 
rived, as described in the Appendix. In Celotti & Ghisellini (2008) 
€b and e c have been found (through spectral modelling) for a sam- 
ple of BL Lacs and FSRQs detected in 7-rays. It was found that 
there exists a good correlation between e c and Lj, and a much more 
scattered correlation (or trend) between es and Lj. Using the rela- 
tions found in that paper, we here set: 



0.1 L,: 



j,45 



e B ~ 0.05 L 



J, 45 



(17) 
(18) 



2.8 Energy of the electrons emitting at the peaks of the SED 

Let us call 7 poa k the random Lorentz factor of the electrons emit- 
ting at the two peaks of the SED. We take advantage of the ob- 
served correlation between 7 poa k and the (comoving) energy den- 
sity l7'(7 P oak) that is the sum of the magnetic and radiation en- 
ergy density integrated up to v — m c c 2 /(/i7 pea k). In other words, 
C^'(7peak) is the sum of the magnetic and the radiation energy den- 
sity available for scattering in the Thomson regime (i.e. efficient 
cooling). This correlation (updated in Celotti & Ghisellini 2008) is 
of the form 



7pcak 



7pcak 



= 10 J 



0.3 



-1/2 



U' ^ 0.3 erg cm 



U' ^ 0.3 erg cm 



(19) 



In Ghisellini et al. (2002) we have interpreted this correlation 
(found by modelling SEDs of blue BL Lac objects) as the result 
of the radiative cooling occurring in one light crossing time. In this 
time electrons cool down to an energy -y c given by 

3m e c 2 



7c = 



oc U" 



(20) 



4<7 T Ai?'£/' 

This scaling stops when 7 C < 71, because in this case the peak 
(in vF v ) is produced by electrons at 71. According to Celotti & 
Ghisellini (2008), this occurs when 7 C ~ 10 3 . For larger energy 
densities, the scaling 7p ea k oc U'~ 1/2 suggests that it is the cooling 
rate (at 7 poa k) which is constant, i.e. 7 oc 7p 0a k^' is tne same for 
all powerful blazars. Note, however, that the scatter in the region of 
large U' is large, making this relation true only approximately. 

2.9 The emitting particle distribution 

The energy distribution of the emitting particles, ^(7) [cm -3 ], is 
assumed to be the same as the one described in Ghisellini, Celotti 
& Costamante (2002) and Celotti & Ghisellini (2008). This distri- 
bution approximates the case of an injection of particles, through- 
out the source, lasting for a finite injection time, equal to t- ln j — 
AR' /c — tpRdiss/c. This is because blazars are variable (flaring) 
sources, and a reasonably good representation of the observed spec- 
trum can be obtained by considering the particle distribution at the 
end of the injection. When the injection stops, particles above 7 C 
have cooled, modifying the energy distribution of the injected par- 
ticles. The latter is assumed to be a broken power-law with slopes 
oc 7 _1 and oc 7~ s below and above the break at 7i n j. We as- 
sume that 2 < s < 3. In the case of fast cooling, occurring when 
7c < 7inj, we have: 



7V( 7 ) 
iV( 7 ) 



oc 7 
oc 7~ 



-(H-i). 



7 > 7inj 

7c < 7 < 7inj 



© 0000 RAS, MNRAS 000, 000-000 



6 G. Ghisellini and F. Tavecchio 



7 < 7c 



(21) 



Since s > 2, in this case the peak energy 7 pca k = 7mj. When 
instead 7 C > 7mj (slow cooling), we have: 

N(j) <x 7" (s+1) ; 7>7c 

AT (7) oc 7 _s ; 7inj < 7 < 7c 

AT (7) tx 7" 1 ; 7 < Tinj. (22) 

For the assumed range of values of s, the peak energy in this case 

is 7pcak = 7c- 

2.10 Summary of free parameters 

Our scheme needs the following parameters: 

• the black hole mass M; 

• the accretion rate Mm, or, equivalently, the disk luminosity in 
units of the Eddington one, Ldisk /-^Edd! 

• the bulk Lorentz factor V; 

• the initial separations of the colliding shells (in the inter- 
nal shock case) or, more generally, the distance, in units of 
Schwarzschild radii, of the dissipation region; 

• the "equipartition" parameters e c and es 

• the jet efficiency factor r7j ct ; 

• the slope s of the injected particle distribution. 

The viewing angle 6 is not considered as a free parameter, since 
we always assume r ~ S and then 9 ~ 1/r. We have rather good 
observational constraints on several of these 9 parameters. For in- 
stance, 10 < r < 20; -Rdiss is a few hundreds Schwarzschild radii; 
the (average) slope of the injected particle distribution is s ~ 2.4- 
2.5; the jet is slightly more powerful than the disk luminosity, lead- 
ing to rij c t ~ 0.2-0.5; as mentioned, we have derived es and e e for 
a number of blazars, and we can use these numbers in our scenario. 
In general, the most important parameters are thus M and Mj n , that 
are the leading quantities characterising our proposed scheme. 



3 SIMPLE CONSEQUENCES 

In this section we derive in an heuristic way a few simple conse- 
quences of our assumptions, leaving for the following sections a 
more detailed description, which requires more technical details. 

• BL Lac/FSRQ division — The first consideration is, in itself, 
one of our assumptions: the BLR exists only when Z/disk/^Edd is 
greater than some value, which we take equal to 3 x 10 -3 . This im- 
mediately implies a division between FSRQs and BL Lac objects, 
defined as objects with and without broad emission lines, respec- 
tively. 

• Existence of "blue" FSRQs — The BLR radius i? B LR oc 
^diak ^ M^ 2 , while the dissipation distance is -Rdiss oc M. 
Therefore there is the possibility that relatively high mass objects, 
with relatively faint disk (but with Ldisk > L c ) have jets pref- 
erentially dissipating beyond the BLR (see also Georganopoulos, 
Kirk & Mastichiadis 2001, Pian et al. 2006). The emitting electrons 
would suffer less radiative cooling, implying a large 7 pca k and a 
"blue" SED. From Eq.[9]and Eq.|4]we have that Rdiss > Rblr for 



0.39M 9 < Ldisk,45 < 0.81M 9 2 



a(T) 



300 



(23) 
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Figure 2. The disk luminosity vs black hole mass plane. Different blazar 
classes occupy different regions of this plane. The dashed line corresponds 
to a disk luminosity equal to a fraction of the Eddington one (for this figure 
we have used 3 X 10 — 3 ). Below this line we assume that the broad line 
region does not exist, or it is very weak. Therefore BL Lac objects live 
below this line: "classic" BL Lac objects just below, and TeV BL Lacs at 
smaller still disk luminosities. According to our scenario, there is another 
important dividing line, corresponding to where i?diss > Ablh [f° r this 
figure we have used a(T) = 300]. Objects having a BLR, but whose jets 
preferentially dissipates beyond it, are "blue" quasars. If blazars can have 
relatively small black hole masses (i.e. even smaller than 10 s Mq), then 
there should be a population of relatively low power quasars, with a red 
spectrum. 



• Existence of "red" FSRQs at relatively low power — This 
is a consequence of considering relatively small black hole masses. 
If FSRQs with -say- M ~ 3 X 10 7 M exist, then they have disk 
luminosities down to Ldisk ~ 1.2 x 10 43 erg s _1 ; and similar jet 
powers (see Fig. 2). Dissipation takes place well within the (small) 
BLR, in regions of large U' implying small 7 pca k and hence a red 
SED. The emitting regions would also have relatively large mag- 
netic fields (Ub oc 1/M 2 , see Eg. I14t. implying large Ub- The 
synchrotron self-Compton emission would then be relatively more 
important than the External Compton component. 

• The Compton dominance — It is defined as the ratio of 
the inverse Compton to synchrotron luminosity Lc/L s and cor- 
responds to (for scattering in the Thomson regime) 



Lc UL + ULt 
L B ~ U B 



(24) 



This defines a "triangular" region in the L^^-M plane, as illus- 
trated in Fig. [2] 



where the external radiation energy density U' c ^ t (in the comoving 
frame) is produced by the BLR and by the dusty torus, while U' s 
is the synchrotron radiation energy density. In general, we expect 
that the Compton dominance is larger for FSRQs with Rdiss < 
Rblr. But since Ub is larger for smaller -Rdiss (smaller masses), 
while ?7c Xt is constant (for equal F), we expect that the Compton 
dominance is on average larger for larger masses. 

In Fig. [2] we illustrate these simple points in the plane Ldisk- 
M, The line dividing BL Lac objects from FSRQ is assumed to to 
be at Ldisk/iEdd = 3 x 10 -3 , the radius of the broad line region 
corresponds to Eq. [4] and the dissipation radius assumes a(T) = 



© 0000 RAS, MNRAS 000, 000-000 



The blazar sequence: a new perspective 7 



300. These values are indicative, but not certain, and there can be a 
substantial scatter. We recall that "red" and "blue" here mean that 
the synchrotron and inverse Compton peak frequencies have low 
or large values, respectively, but not that they are necessarily more 
or less Compton dominated. To make further predictions about the 
predicted SED as a function of M and M- m , we must evaluate the 
Compton dominance in detail. This is done below. 



4 ESTIMATING THE COMPTON DOMINANCE 

In this section we show how it is possible, with a few reasonable as- 
sumptions, to derive the ratio between the inverse Compton to syn- 
chrotron luminosity. This parameter, called Compton dominance, 
will determine in what region of the M -M\ n plane we can find 7- 
ray bright sources (and viceversa where are the 7-ray dim blazars). 

In the Appendix we show the importance of the so-called 
Comptonization y parameter, and how it can be calculated. We de- 
fine it as 



y = axAR' I 7 2 A r (7)d7 



(25) 



In the Appendix we also show that the comoving synchrotron 
and self-Compton radiation energy densities can be expressed as 
Us = yUs and Ug SC — V 2 Ub, while the external Compton en- 
ergy density is U' EC = yU' c ^ t . The asymptotic values of y are: 



1 + 



1 

2 

£c£A U B 



4e c £A 



CD 



1/2 



TP 



17b 



< 1 



> 1 



(26) 



where we have assumed that the scattering process occurs in the 
Thomson regime. In this regime the ratio between the inverse 
Compton and the synchrotron luminosities can be approximated as 



U' 

syn 



UL 



cxt 

Ub~ 



y + 



TP 



(27) 



The first term on the RHS controls the relative importance of the 
SSC component, while the second measures the relative importance 
of the EC component. They are not independent, because the value 
of y is in turn controlled by the (U^/Ub) ratio, as shown by Eq. 
[26] 

Isolating the EC dominance we can focus on FSRQs with 

T^diss < 7?blr.: 



Lbc 



UL 



M 2 
Mm 



«. r 4 [o(r)f 

U B 



M 



Li/Le 



Kr)] 2 

£B 



(28) 



where we made use of Eq. [6] Eq. [9] and Eq. [14] This implies that 
FSRQs of the same mass and increasing M ln form a sequence of in- 
creasing power and decreasing Compton dominance. On the other 
hand FSRQs of increasing mass and constant L, / L-^dd form a se- 
quence of increasing power and increasing Compton dominance. 
Note also the strong dependence on the bulk Lorentz factor. FS- 
RQs of equal masses and accretion rates can form a sequence of 
increasing power and increasing Compton dominance by having 
slightly different values of F. 

The SSC over synchrotron power ratio is simply y. When the 
EC emission is unimportant, y ~ (e c eA/eB) 1//2 , greater than unity 
when e c > e B /tA (Eq.l26l. 
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Figure 3. The disk accretion luminosity (top panel) and the observed bolo- 
metric jet emission (bottom panel), as a function of black hole mass. The 
bulk Lorentz factor is assumed to be 15 for all cases, while we have as- 
sumed a(r) = 300 for FSRQs and 100 for BL Lacs. The solid lines on 
bot panels correspond to L/L Edd = 3 X 10 -3 and on to the region where 
^diss > RbIjH- The different grey levels correspond to different "Comp- 
ton dominance", i.e. different Lq / L E ratios. 



5 RESULTS 

The results of our scheme, for the set of chosen parameters, are 
shown in Fig. [3] and Fig. [4] In this figures the different grey levels 
corresponds to different Compton dominance values, as labelled. 
Consider first Fig. [3] The top panel is equivalent to Fig. [2] Below 
the diagonal line Ldi s k/7/Edd = 3 x 10 -3 we have BL Lacs, above 
the line we have FSRQs. 

The solid lines divide the plane in three regions: in region I, 
above the solid line Ldisk/7/Edd = 3 x 10~ 3 , we have the FSRQs. 
Below this line, in region II, we have BL Lac objects. In region III, 
limited by the two solid lines, we have the "blue" FSRQs. 

The Compton dominance in region I increases for increasing 
mass and decreases for increasing Ldisk for objects with the same 
mass, as described by Eq. 1281 If or a(T) constant]. In region II the 
Compton dominance, determined by y ~ (ectA/eB) 1 ^ 2 , is varying 
much less and is of order unity in the entire region. In region III 
we have the same trend as in region I, but with a lower Compton 
dominance value, determined partly by SSC emission (scaling with 
y) and partly by the torus external IR photons. 

The bottom panel shows the observed bolometric luminosity 
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Figure 4. The peak synchrotron frequency ^ SjPca k (top) and the peak 
Compton frequency v c -p C3 ^ (bottom) as a function of the observed jet lumi- 
nosity. As in Fig. [3] different grey levels indicate different Compton domi- 
nance values. 



produced by the jet as a function of the mass. In this plane objects 
belonging to regions II and III partially overlap, since we can have 
"blue" FSRQs slightly less luminous than BL Lacs. This effect is 
due to the fact that in this illustrative figure we have assumed that 
the dissipation region of BL Lac objects has a(T) = 100, while 
we use a(F) = 300 for FSRQs. Some BL Lac objects, therefore, 
although lacking the BLR, are slightly more efficient (larger ea) 
than "blue" FSRQs, simply because they are more compact. This 
also means that their magnetic field is relatively larger, resulting in 
a less Compton dominated source. 

In Fig. [4] we show the predicted peak frequency of the syn- 
chrotron flux (top panel) and of the IC emission (bottom). Again, 
different grey levels correspond to different Compton dominance 
values. Note that the range of ^ s , pca k and i^ c ,poak values is much 
broader for BL Lac objects than for FSRQs. This is the conse- 
quence of two facts: first, the dependence of 7 pea k on U' is stronger 
for BL Lacs than for FSRQs (7 pca k oc U'" 1 vs 7 pca k oc t/ /_1/ ' 2 
respectively); secondly, the energy density of the external seed pho- 

1/2 

tons is constant, as implied by assuming -Rblr oc L disk and a con- 
stant V: thus when the external radiation energy density dominates, 
it implies a quasi-constant 7 pca k, that in turn makes i/ s , pca k and 
^c, pca k to vary in a small range. 

Fig-E]shows that the maximum values of f s , pea k and I'c.peak 
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Figure 5. Predicted SED as a function of black hole mass and disk lu- 
minosities. For FSRQ (top three thick lines) we have assumed a constant 
L/L Edd ratio (equal to to 0.3) and different black hole masses. For BL 
Lacs, instead (lower four lines thin lines), we the shown SEDs correspond 
to a black hole mass of 10 9 Mq and different L/Z/Edd values, as labelled. 



of FSRQs have a sharp boundary. This is due to the following. 
When the BLR radiation energy density dominates the cooling, we 
have 

3.6 x 10 6 r_B 2.7 x 10 13 rB 

0.944r 2 + B 2 



, s , pcak = 3.6X10 rB 7pcak = (17/12)[/BLRr2 + [/s 

where we have used Eq. [6] and Eq. [19] The maximum ^ s , pca k is 
reached for B 2 = 0.944r 2 . In our case, since T = 15, we have 
fs.peak = 1.4 X 10 13 Hz. 

Blue FSRQs, on the other hand, suffer much less cooling 
(they dissipate beyond the BLR) and consequently their 7 pea k and 
z^ijeak are larger, explaining the existence of the gap (in i-'s.pea.k) 
of the two populations. 

It is worth to stress that blazars will not all have the same value 
of r (that can change even in a single source). The sharp boundaries 
described above and shown in Fig.|4]will disappear once we allow 
for a distribution of F-values. 



5.1 Predicted SED 

In Fig. [5] we show several predicted SEDs of FSRQs and BL Lac 
objects. For FSRQs the spectra correspond to objects with the same 
ratio Z/disk/^Edd = 0.3 and different masses, while for BL Lacs 
we have the same mass (10 9 Mq) and different Z/disk/^Edd ra- 
tios. In agreement with Fig. [3] FSRQs are more Compton domi- 
nated as the mass (and then the bolometric luminosity, for constant 
-f'disk/i'Edd) increases, with almost constant peak frequencies. 
Note that, increasing the mass, the IC luminosity increases, while 
the synchrotron component decreases. This can be understood re- 
calling that L' s oc yUBRdi ss and inserting the appropriate scaling 
of Ub (Eq. [14) and y (Eq. [35j, yielding L' s oc e c esL 2 /i? diss oc 



Hz(2S 
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Figure 6. Predicted SED as a function of the bulk Lorentz factor for a 
fixed black hole mass (1O 9 M0) and for L/L^dd — 0.1 (top panel) and 
L/Lsdd = 10 — 2 (bottom panel). For the latter case, note the SED of the 
r = 16 case, remarkably different from the others. In this case i?diss > 
i?BLR an d no external BLR photons contribute to the radiative cooling. 
There is instead a contribution coming from the torus. The reduced cooling 
implies a large 7 pca k, making the SED similar to the one of a blue BL Lacs. 
However, the broad lines do exists, and the objects would be classified as 
"blue" quasars (with a synchrotron peak frequency ^ pca k ~ 10 15 Hz). 



e c es(Ldisk/-£'Edd) 2 - Therefore, for a constant ratio Ldisk/^Edd, 
the synchrotron luminosity L s oc e c es which we assume to scale 



as L . 



-(0.4+0.25) 



(see Eq.[T7]and Eq.gg). 



The spectral indices in the ~ GeV (a 7 ) and X-ray (ax) bands 
anti-correlate: at low luminosities we have a steep ax (due to the 
tail of the synchrotron emission) and a hard a 7 (due to a still rising 
SSC spectrum), while at large powers we have a flat ax (due to 
the low energy part of the EC spectrum) and a steep q 7 (due to the 
high energy tail of the same component). This is in agreement to 
what observed (e.g. Comastri et al. 1997). 

The series of SED in Fig.[5]resembles the phenomenological 
sequence of Fossati et al. (1998). In order to reproduce the observed 
increase of the Compton dominance with the total power for FSRQs 
we assume the same Ldisk/^Edd and change the black hole mass. 
This choice is dictated by the fact that, in the present scheme, the 
Compton dominance increases with the mass (as clearly visible in 
Fig.fJJ. Fixing the mass and increasing the disk luminosity, instead, 
would have the effect to reduce the Compton dominance, contrary 
to the observed trend. Low power BL Lac objects, instead, display 
in the blazar sequence a more or less constant Compton dominance 
and we can reproduce this branch of the blazar sequence using a 
fixed (large) black hole mass and varying the accretion rate. 



6 CHANGING PARAMETERS 

In the previous section we have shown the SED expected for differ- 
ent values of the mass and accretion rate. Here we briefly discuss 
the consequences of changing some other input parameters of the 
model. 

• Changing the bulk Lorentz factor F. In Fig.|6]we show the pre- 
dicted SED as a function of the bulk Lorentz factor for L/Z/Edd = 
0.1 (top panel) and L/Lsdd = 10~ 2 (bottom panel). The black 
hole mass is kept constant (M = 10 9 Mq). The bulk Lorentz fac- 
tors are in the range 7-16. It can be seen that even small changes of 
r imply large changes in the observed SEDs, as described by Eq. 
1281 as long as i?di ss < J^blr- Since we have assumed a(F) oc T 2 
for this test, we have L C /L B oc T 8 . The top panel shows a mono- 
tonic sequence, since for all T we have i?di ss < -Rblr- The bottom 
panel, instead, shows the case of a "blue" quasar: for the largest 
adopted T, in fact, we have i?di ss > Rblh., and the external pho- 
tons contributing to the EC spectrum are coming from the torus 
only, contributing to the comoving radiation energy density accord- 
ing to Eq.[8] This implies a slower electron cooling, a larger 7 pca k 
and consequently larger peak frequencies. Observers would how- 
ever detect the broad emission lines in this blazars, classifying it as 
a "blue" FSRQ, with a synchrotron peak frequency f pca k ~ 10 15 
Hz. 

• Changing e c and e_g. There are two different regimes: when 
the external seed photons dominate, e c controls the total produced 
luminosity, but not the Compton dominance, which is instead pro- 
portional to e^ 1 (see Eq|28ll. When instead the high energy spec- 
trum is dominated by the SSC flux, the Compton dominance is pro- 
portional to (ee/eB) 1 ^ 2 (as in Gamma Ray Burst afterglows). 

• Changing i?diss- Although we always assume Rdi ss oc M, we 
can change the constant of proportionality [a(r)]. As a rule, if we 
are inside the BLR, a smaller iidiss means a larger magnetic field, 
and thus a smaller Compton dominance (see Eq.[l4] Eq.|28t, Also, 
increasing a(T), makes the case i?di ss > Rblr. more probable, 
enlarging the parameter space for blue quasars. 

• Changing s. Besides changing the predicted slope above the 
peaks, it changes e\ for objects emitting in slow cooling: the 
steeper s, the less efficient the source (small ea). Since the major- 
ity of blazars emitting in the slow cooling regime are the ones with 
no BLR (and torus), a steeper s also means less (self-) Compton 
dominance. 

• Changing r)j . This parameter act as a normalisation: larger ryj 
means that the jet is more powerful overall. Keeping all other pa- 
rameter constant, increasing r)j increases the value of the magnetic 
field, and therefore it decreases the Compton dominance when the 
external seed photons dominate the radiation energy density. When 
they do not, changing rjj has no effect on the Compton dominance. 



7 CAVEATS 

Our proposed scenario is necessarily highly idealised and is aimed 
at describing the averaged behaviours of different classes of blazars 
and also of a single blazars that can have different states. In this 
section we therefore try to summarise a few caveats that should be 
kept in mind. 

• Leptonic model. Probably the most basic assumptions of ours 
is that the emission comes from the acceleration and the radiative 
cooling of leptons, and we have neglected the possible presence of 
relativistic hadrons (Mannheim 1993, Aharonian 2000, Miicke et al 
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2003). While the variability patterns observed in some (admittedly, 
a few) blazars seem to favour the leptonic model, the hadronic sce- 
nario is not ruled out yet. 

• Single zone emission region. We have assumed that most of 
the SED is produced in a single region, characterised by a single 
valued and uniform magnetic field, with a single particle distribu- 
tion, and so on. This is the most commonly adopted simplification, 
that received support from the few cases in which we could see co- 
ordinated variability in different energy bands, but consider that i) 
there surely are several emitting regions, contributing differently to 
the entire SED. The region producing the radio flux, for instance, 
is surely much larger than what we consider here, and could con- 
tribute (by its SSC flux) to the X-ray radiation; ii) the jet could be 
structured, composed by a fast spine surrounded by a slower sheat- 
layer, as envisaged in Ghisellini, Tavecchio & Chiaberge (2005). 
Some support in this sense for TeV BL Lac objects comes from the 
fact that these sources have "slow" VLBI knots, often moving sub- 
luminally (Piner & Edwards 2004, Piner, Pant & Edwards 2008 
and references therein), and from direct observations of some radio 
edge-brightening (Giroletti et al. 2004, 2006). 

• A different i?di ss for BL Lacs. Associated to the previous point 
is the possibility that not all blazars have the same a(T) appearing 
in the relation R = a(T)R s . We have in fact assigned, to BL Lac 
objects, a different a(T) than in FSRQs, but the real value could be 
different from what we have adopted. 

• Different bulk Lorentz factors F for different classes of ob- 
jects. There is the need, for TeV BL Lacs, of a very large T, larger 
than for other blazars. The reasons are the very short variability 
timescales observed for the TeV flux (3 minutes in PKS 2155-304, 
Aharonian et al. 2007; see the implications in Begelman, Fabian & 
Rees 2008 and models in Ghisellini & Tavecchio 2008); and the 
necessity, in single zone SSC models, to account for the large sep- 
aration, in frequency, of the two peaks of the SED (e.g. Konopelko 
et al. 2003; Tavecchio, Maraschi & Ghisellini 1998). 

• The relation y pe ak(U). This seems well obeyed by all ob- 
served sources of low-medium power, while for powerful objects 
we have a larger scatter. There is a branch at large U for which 
7 poak oc U^ 1 instead of cx U~ 1/2 (see Celotti & Ghisellini 2008). 

• Thomson vs Klein-Nishina scattering regime. Our estimates 
assume that the inverse Compton process occurs in the Thomson 
regime, so only for 7 pea k not extremely large. We can estimate how 
large by the following argument. If the EC process is important, 
then cooling is efficient, 7 pca k < 10 3 -10 4 and the bulk of scatter- 
ings indeed occurs in the Thomson regime. When instead the SSC 
process is the main producer of high energy photons, then we may 
neglect Klein Nishina effects when 7 pea k lw' a /(m e c 2 ) < 1, where 

~ 3.6 x 10 6 i?7 pcak is the synchrotron peak frequency in the co- 
moving frame. Therefore the Klein-Nishina regime is unimportant 
if 

>-<(3^?i^) 1/3 - 3 - 2x 1045-1/3 (30) 

This corresponds to observed synchrotron peak frequencies v B ~ 
3.7 x 10 15 B 1 ' 3 S Hz and self Compton frequencies u c = 
Tpeak^s ~ 3.8 x 10 24 B- 1/3 S Hz ~ 15.6B~ 1/3 (5 GeV With 
5 ~ 10 and B ~ 0.1 G, blazars with Compton peak energies 
below ~ 70 GeV should still be in the Thomson limit. 

• SED based on flares. Our scenario is inevitably built upon 
what we already know, and this is mostly limited by the sensitivity 
of EGRET. Therefore, the complete SED we could construct, from 
the radio to the GeV band, are most likely associated to sources in 



very active (7-ray) states, and might not be indicative of the more 
standard, but more persistent, blazar spectrum. 

• Viewing angle. The models presented here assume that the 
jet is observed with a small angle 9 ~ 1/r ~ 4 deg. Accord- 
ing to the current unification schemes, besides these highly beamed 
sources a large number of moderately less beamed blazars should 
exist (9 ~4-7 deg), for which the non-thermal continuum of the 
jet is less amplified and thus less luminous. An effective way to 
distinguish between these sources and true low power FSRQs is 
through the level of the optical emission, that in the misaligned ob- 
jects should be dominated by the "naked" thermal disk emission 
(see also Landt et al. 2008). 



8 IMPLICATIONS FOR COSMIC EVOLUTION 

The association of different classes of blazars with the mass of the 
black hole and with the accretion rate can have important effects on 
the predicted evolution of blazars. 

Black hole of large masses (M > 1O 9 M ) are rare, while 
black holes of smaller masses are more common (at any redshift). 
If the accretion rate in blazars evolves in redshift, as indicated by 
the redshift evolution of radio-quiet quasars, we should have larger 
accretion rates in the past. As a consequence, blazars whose accre- 
tion rate is less than a critical value - sources that we associate with 
BL Lac objects - should become rarer in the past, implying a neg- 
ative evolution. On the contrary, FSRQs should be more common 
in the past, since a larger fraction of massive black holes accretes 
at rates larger than the critical one. 

At large redshifts, when the black hole is fully accreting, we 
have powerful disk and jets and a BLR located at relatively large 
distances. These are the conditions to have a powerful, Comp- 
ton dominated, FSRQ. The same object will lower its accretion 
rate with time, the disk becomes a less efficient radiator, the BLR 
shrinks, and there should be the possibility to have a transitional 
object, with emission lines of moderate luminosity, and with a 
SED that can be produced sometimes within, sometimes outside, 
the BLR region (i.e. the red/blue FSRQ case in Fig. [2}. When the 
accretion rate decreases below the critical value (i.e. Li/I/Edd < 
3 x 10~ 3 , the accretion changes mode, becoming radiatively ineffi- 
cient. The corresponding BLR (if it exists at all) becomes weak and 
very close to the black hole. Jet dissipation occurs always beyond 
the BLR, and the source becomes a BL Lac object. Synchrotron 
and SSC are the main emission processes, and the position of the 
corresponding peaks in the SED depends on the amount of radiative 
cooling (bluer for less cooling). 

These simple ideas are in line with the scenario proposed by 
Cavaliere & D'Elia (2002), in which blazars are born initially pow- 
erful (high accretion rate), and transform themselves into BL Lac 
objects when the accretion rate decreases. The evolution of blazars 
as a whole should be the same of radio-quiet quasars (Maraschi & 
Rovetti 1994), but with black hole masses larger than ~ 3xl0 7 M o 
(if this is the real limit to have a radio loud object). SED-wise, the 
evolution is more complicated, since for low masses and interme- 
diate accretion rate, a FSRQ can be blue (i?di ss > Rblr,), but still 
have visible broad lines. 



9 CONCLUSIONS 

We have proposed an extension of the scheme based on the so- 
called "blazar sequence", taking advantage from the new facts and 
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ideas discovered or suggested since the "old" blazar sequence was 
published. Our aim is not to propose a new model or scenario for 
blazars, but rather to view them from a different point of view, 
enabling a different perspective. Technically, we pass from a one- 
parameter (observed bolometric luminosity) to a two-parameter se- 
quence (M and M; n ), but, more physically, we try to link the prop- 
erties of the jet to the main properties of the accretion disk. To do 
this, we had to make a number of simplifying assumptions, that we 
think are rather reasonable or even supported by strong observa- 
tional evidence, but of course there remains some uncertainties. In 
addition to this, we stress that, at this stage, we can only (and nec- 
essarily) describe averaged properties of blazars and some trends. 
But the link between the two fundamental parameters (M and M ln _, 
or, rather, Ldisk and M) characterising the new sequence and the 
blazar SED makes our proposal perfectly testable (and falsifiable). 
The main results are: 

• Blue quasars may exist. If the dissipation region of the jet is 
beyond the broad line region, the radiation energy density seen in 
the comoving frame is not large, resulting in less electron cooling, 
and thus large electron energies. This in turn results in large syn- 
chrotron and inverse Compton peak frequencies. This might hap- 
pen for large black hole masses and intermediate accretion rates 
(for small M there are no broad lines anymore, and the objects be- 
comes a BL Lac, for large M the radius of the BLR is large, making 
it difficult that the jet produces most of its luminosity beyond it). 
If the dissipation region is related to the value of the bulk Lorentz 
factor, as in the internal shock scenario, then an additional require- 
ment is that T should be relatively large. Note that the region of the 
parameter space corresponding to blue quasars is not large, and its 
extension depends upon V and the relation between the radius of 
the BLR and the disk luminosity. At this stage, it is then difficult to 
quantify more precisely how many blue emission line blazars there 
are. 

• Low power blazars (red FSRQs; red and blue BL Lacs) should 
exist. They correspond to blazars associated with low mass black 
holes. The blue BL Lacs, corresponding to smaller accretion rates 
and hence smaller jet powers, can be missed more easily (than red 
blazars) in flux limited samples. It is possible that this effect can 
account for the large number of red blazar in deep flux limited sam- 
ples (see e.g. Padovani et al. 2007) especially if blazars with black 
holes of small masses are more numerous. To demonstrate this, we 
of course need detailed simulations, that we leave for future work. 

• Red low power blazars should exist. Again they simply cor- 
respond to objects with black holes of relatively small masses and 
relatively large accretion rates, hence with a "normal" broad line 
region. Low power red sources should be easily distinguished from 
the powerful misaligned counterparts (expected to have a similar 
SED) since in these latter objects the powerful disk emission is ex- 
pected to dominate over the non-thermal continuum in the optical- 
UV region (see Landt et al. 2008). 

• Not only the power of jets, but also the SED that they produce 
is linked with M and Mi n . This is the most important result of our 
study. 

To confirm or falsify our ideas we need to know the SED, the 
mass of the black hole, and the disk luminosity (or, equivalently, 
the luminosity of the broad emission lines). We hope to soon enter 
in the "GLAST era", greatly helping to determine the blazar bolo- 
metric luminosity and their high energy SED. If the blazar sample 
of EGRET counts ~60 objects, a improvement of a factor 10 in 
sensitivity means roughly that GLAST will detect ~ 2000 blazars. 
Assuming that Swift can observe, simultaneously and with its op- 



tical and X-ray telescopes, 1/10 of them, we can have a sample of 
~200 blazars with an optical to GeV simultaneous coverage. This 
means a large improvement in number, and a huge improvement in 
quality, of our knowledge. These blazars can form a "golden" sam- 
ple for ground based follow up, aiming to measure (if there are) the 
luminosity the broad emission lines and estimate the mass of the 
central black hole. With these information, we can easily confirm 
or reject the new blazar sequence we propose. 
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Appendix 

The Comptonization y parameter 

Our method uses the fact that we "know" the total radiation 
energy density produced in the dissipation region, since it 
corresponds to the fraction of energy lost by the electrons in 
one light crossing time: TJ[ = tpJJ c . This radiation energy 
density U' T is the sum of the synchrotron and inverse Comp- 
ton emission^. 

Consider first the synchrotron emission, with 
a corresponding cooling rate given by j s ~ 
(4/3)crTcJ7s7 2 /(m c c 2 ). Neglecting self-absorption, 
the total (comoving) synchrotron luminosity is 



L' s = Trr 2 AR'm c c 2 / iV (7)7^7 



(31) 



We now divide and multiply by n = / AT (7)^7 and define 
r = nAR'n. We can also write L' as Ulnr 2 c. We then have 



K = u B gr(7 2 ) = yU E 



(32) 



where the last equivalence corresponds to define the Comp- 
tonization parameter y = (4/3)t(j 2 ). Assuming that the 
scattering process occurs in the Thomson regime we can 
analougously set [7£ sc = yU' s = y 2 U B , and U^ c = yU' cyit . 
The sum of these components constitutes (7 r ': 

U' t = U' syn + U^ SC + U EC = yU B +y 2 U B + yU^ t (33) 

Dividing both parts by Ub we obtain an equation dependent 
on the microphysical parameters es, e c and ea (defined in 
§2.7, see Eqs. [T4l|T5l >, and dependent on Lj only through the 
ratio U' cy . t /U B : 



EL 
U B 



£A£c 



it + y 1 + 



UL 



U B 



The solution is: 
1 

V = 2 



TP 
U B 



4e c £A 
£B 



1/2 



(34) 



I + S* ! > (35 i 



For f/gxt/^B much smaller or greater than unity we have the 
limiting cases of Eq. [26] 

Finding ea 

To find ea, consider an injected relativistic electron distribu- 
tion that has a power law form, (5(7) = Qo7 s > between 
some 71 and 72. The fraction 6a of radiated energy (in one 
light crossing time) can be approximated as 



£A = 



P 2 -y 1 

J 7pe a k ' 



5 ^7 / 7l 



S^'T 1 -'dry V7 P oak 



s-2 



s > 2 



(36) 



The value of 7 poa k depends upon U' = Ub +' s +U^t (Eq- 
n~9b . In particular, the synchrotron radiation energy density 



3 U' T should not be confused with U' of Eq.fl9] which is the magnetic plus 
radiative energy density available for efficient scattering and contributing to 
the cooling. 
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depends upon y (Eq.[33l>. i.e. it depends on ea itself. Despite 
this circularity problem, we can solve the system iteratively, 
until convergence (obtained in a few cycles). 
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